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We consider the magnetorotational instability (MRI) of a hydrodynamically stable Taylor-Couette 
flow with a helical external magnetic fleld in the inductionless approximation deflned by a zero 
magnetic Prandtl number (Pm = 0). This leads to a considerable simplification of the problem 
eventually containing only hydrodynamic variables. First, we point out that the energy of any 
perturbation growing in the presence of magnetic field has to grow faster without the field. This is 
a paradox because the base flow is stable without the magnetic while it is unstable in the presence 
of a helical magnetic fleld without being modifled by the latter as it has been found recently by 
HoUerbach and Riidiger [Phys. Rev. Lett. 95, 124501 (2005)]. We revisit this problem by using a 
Chebyshev collocation method to calculate the eigenvalue spectrum of the linearized problem. In 
this way, we conflrm that MRI with helical magnetic fleld indeed works in the inductionless limit 
where the destabilization effect appears as an effective shift of the Rayleigh line. Second, we integrate 
the linearized equations in time to study the transient behavior of small amplitude perturbations, 
thus showing that the energy arguments are correct as well. However, there is no real contradiction 
between both facts. The linear stability theory predicts the asymptotic development of an arbitrary 
small-amplitude perturbation, while the energy stability theory yields the instant growth rate of 
any particular perturbation, but it does not account for the evolution of this perturbation. Thus, 
although switching off the magnetic fleld instantly increases the energy growth rate, in the same 
time the critical perturbation ceases to be an eigenmode without the magnetic fleld. Consequently, 
this perturbation is transformed with time and so looses its ability to extract energy from the base 
flow necessary for the growth. 
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The magnetorotational instability (MRI) is thought to 
be responsible for the fast formation of stars and entire 
galaxies in accretion disks. For a star to form, the matter 
rotating around it has to slow down by transferring its 
angular momentum outwards. Without MRI this pro- 
cess would take much longer than observed because the 
velocity distribution in the accretion disks seems to be 
hydrodynamically stable while the viscosity alone is not 
sufficient to account for the actual accretion rates. It 
was suggested by Balbus and Hawley Q that a Kep- 
lerian velocity distribution in an accretion disk can be 
destabiHzed by a magnetic field analogously to a hy- 
drodynamically stable cylindrical Taylor-Couette flow as 
it was originally found by Velikhov (21 and later anal- 
ysed in more detail by Chandrasekhar In this case, 
the effect of "frowziness" of the axial magnetic fleld in a 
well conducting fluid provides an additional coupling be- 
tween the meridional and azimuthal flow perturbations 
that, however, requires a magnetic Reynolds number of 
Rm ~ 10. For a liquid metal with the magnetic Prandtl 
number Pm ~ 10~^ — 10^^ this corresponds to a hydro- 
dynamic Reynolds number Re — Rm/Pm ~ 10^ — 10^ 
Thus, this instability is hardly observable in the 
laboratory because any conceivable flow at such Reynolds 
number would be turbulent. However, it was shown re- 
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cently by Hollerbach and Riidiger 01 that MRI can take 
place in the Taylor-Couette flow at Re ~ 10'^ when the 
imposed magnetic field is heHcal. The most surprising 
fact is that this type of MRI persists even in the induc- 
tionless limit of Pm — where the critical Reynolds num- 
ber of the conventional MRI diverges as ~ 1/Pm. This 
limit of Pm = formally corresponds to a poorly con- 
ducting medium where the induced currents are so weak 
that their magnetic field is negligible with respect to the 
imposed field. Thus, on one hand, the imposed magnetic 
field does not affect the base fiow, which is the only source 
of energy for the perturbation growth. But on the other 
hand, perturbations are subject to additional damping 
due to the Ohmic dissipation caused by the induced cur- 
rents. 

We show rigorously that the imposed magnetic field 
indeed reduces the energy growth rate of any particular 
perturbation. On one hand, this impHes that the energy 
of any perturbation, which is growing in the presence 
of magnetic field, has to grow even faster without the 
field and vice versa. But on the other hand, the fiow 
which is found to be unstable in the presence of magnetic 
field is certainly known to be stable without the field. 
This apparent contradiction constitutes the paradox of 
the inductionless MRI which we address in this study. 

Consider an incompressible fiuid of kinematic viscos- 
ity ly and electrical conductivity a filling the gap between 
two infinite concentric cylinders with inner radius Ri and 
outer radius Ro rotating with angular velocities fli and 
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rio, respectively, in the presence of an externally imposed 
steady magnetic field Bq = B^e^ + B^e^ with axial and 
azimuthal components B-, ~ Bq and — PB^Ri/r 
in cylindrical coordinates (r, 0, z), where /3 is a dimen- 
sionless parameter characterizing the geometrical helic- 
ity of the field. Further, we assume the magnetic field 
of the currents induced by the fluid flow to be negligi- 
ble relative to the imposed field that corresponds to the 
so-called inductionless approximation holding for most 
of liquid-metal magnetohydrodynamics characterized by 
small magnetic Reynolds numbers Rm = ^oavgL ^ 1, 
where /io is the magnetic permeability of vacuum, vq and 
L are the characteristic velocity and length scale. The 
velocity of fluid flow v is governed by the Navier-Stokes 
equation with electromagnetic body force 

dv 



dt 
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where the induced current follows from Ohm's law for a 
moving medium 



j = ct(E + v X Bo) 



(2) 



In addition, we assume that the characteristic time of 
velocity variation is much longer than the magnetic dif- 
fusion time To 3> Tm = p-qctL^ that leads to the quasi- 
stationary approximation, according to which V x E = 
and E = — V$, where $ is the electrostatic potential. 
Mass and charge conservation imply V • v = V • j = 0. 

The problem admits a base state with a purely az- 
imuthal velocity distribution vo(r) = e^vo{r), where 
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Note that the magnetic fleld does not affect the base flow 
because it gives rise only to the electrostatic potential 
*i'o('") — Bq J vo{r)dr whose gradient compensates the 
induced electric fleld so that there is no current in the 
base state (jg = 0). However, a current may appear in a 
perturbed state 
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where vi, pi, j^, and $1 present small-amplitude pertur- 
bations for which Eqs. I^Q El after linearization take the 
form 

(9vi 

— — -I- (vi • V)vo + (vo • V)vi 

ot 
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= — Vpi + i^VVi + -ji X Bo (3) 
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In the following, we focus on axisymmetric perturba- 
tions, which are typically much more unstable than non- 
axisymmetric ones 0- In this case, the solenoidity con- 
straints are satisfled by meridional stream functions for 
fluid flow and electric current as 

V = ve^ + V X (tpe^), j = je^ + V x (he^). 



Note that h is the azimuthal component of the induced 
magnetic fleld which is used subsequently as an alterna- 
tive to $ for the description of the induced current. In 
addition, for numerical purposes, we introduce also the 
vorticity a; — ijje^ -I- V x (we^) = V x v as an auxiliary 
variable. Then the perturbation may be sought in the 
normal mode form 



{vi,uJi,^i,hi} (r,t) = |-e,(2;,-0,/i| (r) 
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where 7 is in general a complex growth rate and k is 
the axial wave number. Henceforth, we proceed to di- 
mensionless variables by using Ri, Rl^/v, Ri^i, Bq, and 
aBoRifli as the length, time, velocity, magnetic fleld, 
and current scales, respectively. The nondimensionalized 
governing equations read as 



jv = JJkV + iieik (r llj r tp + tia iKh, [b) 

7W = DkLO + 2Reiknv-'Ra^ik(^ikii; + 2Pr'''^hj ,{6) 

= Dk^ + uj, (7) 

= Dkh + ik (^i) ~ 2l3r-'^ip) , (8) 



where Dkf = r ^ ('"/')' — (r ^ + k^)f and the prime 
stands for d/dr; Re = R^fli/v and Ha — RiBo^/a/ {pv) 
are Reynolds and Hartmann numbers, respectively; 
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is the dimensionless angular velocity of the base flow 
deflned using A — Ro/Ri and p, = flo/^i- Boundary 
conditions for the flow perturbation on the inner and 
outer cylinders at r = 1 and r = A, respectively, are 
{} = — ip' — 0. Boundary conditions for h on insulat- 
ing and perfectly conducting cylinders, respectively, are 
/i = and {rhy = at r = 1; A. 

The governing Eqs. for perturbation amplitudes 

were discretized using a spectral collocation method on 
a Chebyshev-Lobatto grid with a typical number of in- 
ternal points N = 32 — 96. Auxiliary Dirichlet bound- 
ary conditions for uj were introduced and then numer- 
ically eliminated to satisfy the no-slip boundary condi- 
tions ip' — 0. Electric stream function h was expressed in 
terms of v and ip by solving Eq. ^ and then substituted 
in Eqs. lO that eventually resulted in the 2N x 2N 
complex matrix eigenvalue problem which was efflciently 
solved by the LAPACK's ZGEEV routine. 

In addition, Eqs. were discretized by using a 

Chebyshev-tau approximation and integrated forward in 
time by a fully implicit 2nd order scheme with linear ex- 
trapolation of convective and magnetic terms. We tested 
the numerical code by flnding a few leading eigenmodes 
and eigenvalues by the so-called "snapshot" method 
and compared to the results of the above described code 
as well as to the linear instability results and 
Agreement was at least three signiflcant digits. 
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Figure 1: Critical Reynolds number versus /i for insulating 
cylinders with A = 2 at various helicities /3 and fixed Hart- 
mann number Ha = 15. 



Figure 2: Limiting value of versus the helicity /3 for insulat- 
ing and perfectly conducting cylinders with A = 2 at various 
Hartmann numbers. 



Equations (0 © straightforwardly lead to the kinetic 
energy variation rate of a virtual perturbation Vi satis- 
fying the incompressiblity constraint and the boundary 
conditions. Multiplying Eq. ^ scalarly by Vi and then 
integrating over the volume V which extends axially over 
the perturbation wavelength, we obtain 



dEi 
dt 



[(vi • V)vi] -VodV 



a 



where Ei — ^ J 'v\dV is the energy of perturbation. The 
first integral in the equation above accounts for the in- 
teraction of the perturbation with the basic flow which 
is not affected by the magnetic fleld as noted above. The 
sign of this integral may vary depending on vi. Thus, 
this term presents a potential source of energy. In con- 
trast, the second term is negative deflnite presenting an 
energy sink due to both viscosity and conductivity. Since 
the current is induced only in the presence of a magnetic 
field while the source term does not depend on the mag- 
netic field, we conclude that the instant growth rate of 
any given perturbation has to be lower with magnetic 
field than without it 
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The following results concern cylinders with A = 2, 
as in Ref. 0|. As seen in Fig. which shows the 
critical Reynolds number as a function of /i for Hart- 
mann number Ha =15 and various geometrical helici- 
ties f}, the Hnear instability threshold can indeed extend 
well beyond the Rayleigh line jic = = 0.25, defined 
by ci(r^ri) /dr = 0, when the magnetic field is helical 
(/3 7^ 0). In contrast to Pm ^ the range of in- 
stability is Hmited by ^max, which is plotted in Fig. [21 
depending on the geometrical helicity (3 at various Hart- 



mann numbers Ha for both insulating and perfectly con- 
ducing cylinders. The critical Re tends to infinity as /i 
approaches /imax as in the nonmagnetic Taylor-Couette 
instability. Thus, in the inductionless approximation, the 
destabiHzing effect of a helical magnetic field appears as 
a shift of the Rayleigh fine towards higher ^. The shift is 
especially pronounced for perfectly conducting cylinders 
at « 4. 

The results of time-integration of the Hnearized prob- 
lem are illustrated in Fig. for a perturbation with 
/c = 2 at Re = 2000. This perturbation is unstable in the 
presence of a magnetic field with [3 = A and Ha = 15 
(Rcc = 1554, kc = 2.5) and stable without the field 
because ^ = 0.27 > /ic. First, we integrate an arbi- 
trary, sufficiently small initial perturbation for a suffi- 
ciently long time so that the unstable mode dominates 
but still remains small for the linear approximation to be 
valid. Then we "switch off" the magnetic field by setting 
Ha = 0. Note that we assume the field to be instantly 
absent when it is switched off. So we just compare the 
evolution of the given perturbation with and without the 
field. As seen on the first inset of Fig. El the energy 
of an unstable perturbation indeed starts to grow faster 
instantly after the magnetic field is switched off. How- 
ever, the growth keeps only for a short time and then 
the energy quickly decays as predicted by the linear sta- 
bility analysis. Note that the energy keeps decaying in 
an oscillatory way because the dominating perturbation 
without the field is not a pure traveling wave but rather 
a superposition of two oppositely traveling waves which 
both have the same decay rate and frequency whereas 
the amplitude ratio of both waves is determined by the 
initial condition. 

The magnetic field is switched on again at the instant 
t = 0.1. The corresponding evolution of the perturba- 
tion energy is shown on the r.h.s. of Fig. Elin enlarged 
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Figure 3: Time evolution of the energy of the dominating 
perturbation with = 2 at Re = 2000 after switching the 
magnetic field off and later on again for /i — 0.27, Ha = 15, 
and /3 = 4. Extrapolation shows how the evolution would 
proceed without the change of the magnetic field. 

scale. As seen in the second inset, the energy decay rate 
instantly increases in accordance to ijHJ when the mag- 
netic field is switched on. However, after a short transient 
the perturbation energy resumes the growth in agreement 
with the linear stability analysis. In this case, the energy 
growth is purely exponential because the dominating per- 
turbation is a single traveling wave. 

Thus, this particular example of time integration con- 
firms the validity of Eq. ^ which applies in general 
to any arbitrary perturbation. The energy of an unsta- 
ble perturbation indeed starts to grow faster when the 
magnetic field is switched off. However, there is no real 
contradiction with the linear stability predictions because 
the energy grows only for a limited time and then turns 
to decay as predicted by the Hnear stability. It is impor- 
tant to stress that the linear stability theory predicts the 
asymptotic development of an arbitrary small-amplitude 
perturbation, while the energy stability theory yields the 
instant growth rate of any particular perturbation, but it 
does not account for the evolution of this perturbation. 
Thus, although switching off the magnetic field instantly 
increases the energy growth rate of the most unstable 



as well as that of any other perturbation, in the same 
time the critical perturbation ceases to be an eigenmode 
without the magnetic field. Consequently, this pertur- 
bation is transformed with time and so looses its ability 
I to extract energy from the base fiow necessary for the 
growth. Analogously, switching on the magnetic field 
fe causes an instant decrease of the growth rate of any par- 
g ticular perturbation because of Ohmic dissipation, while 
I the magnetic field transforms the perturbation so that it 
I becomes able to extract more energy from the base fiow 
and so eventually grows. 

To understand the physical mechanism of this instabil- 
ity, note that a heHcal magnetic field, in contrast to pure 
axial or azimuthal fields, provides an additional coupling 
between meridional and azimuthal fiow perturbations. In 
a helical magnetic field with axial and azimuthal compo- 
nents, the radial component of the meridional fiow per- 
turbation induces azimuthal and axial current compo- 
nents, respectively. Interaction of this current with the 
imposed magnetic field results in a purely radial electro- 
magnetic force which retards the original perturbation. 
So, it has a stabilizing effect similar to the radial deforma- 
tion of magnetic fiux lines in the conventional MRI 
However, in the perturbation of finite wavelength there 
is also a radial current component associated with the 
axial one as required by the solenoidity constraint. This 
radial current interacting with the axial component of 
the heHcal magnetic field gives rise to the azimuthal elec- 
tromagnetic force, thus coupling the meridional and az- 
imuthal fiow perturbations similarly to the conservation 
of the angular momentum in the purely hydrodynamic 
Taylor-Couette instability or the azimuthal twisting of 
axial magnetic fiux lines in the conventional MRI. Note 
that the latter effect also renders the imposed axial mag- 
netic field locally heHcal that, however, requires Pm > 
and Re 1 /Pm. When the imposed magnetic field is he- 
lical, the inductionless approximation defined by Pm = 
is applicable to MRI where it leads to a considerable sim- 
plification of the problem containing only hydrodynamic 
variables as in the classical Taylor-Couette problem. 
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